Possible heavy molecular states composed of a pair of excited charm-strange mesons 
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The P-wave charm-strange mesons D a o (2317) and D a \ (2460) lie below the DK and D* K threshold 
respectively. They are extremely narrow because their strong decays violate the isospin symmetry. 
We study the possible heavy molecular states composed of a pair of excited charm strange mesons. 
As a byproduct, we also present the numerical results for the bottonium-like analogue. 
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I. INTRODUCTION 



In the past seven years, the experimentally observed new charmonium or charmonium-like states include A (3872), 
y(3940), F(4260), Z(3930), A(3940), Y(4325), F(4360), F(4660), Z+(4430), Z+(4050), Z+(4250) and F(4140) etc 
flll| - It's difficult to accommodate all these states especially those charged ones in the conventional quark model. 
II • Many of these new states lie close to the threshold of two charmed mesons. A natural speculation is that some of 
them may be the molecular states composed of two charmed mesons 

In the framework of the meson exchange model, we have investigated the possible loosely bound molecular states 
composed of a pair of the ground state S-wave heavy mesons and a pair of S-wave and P-wave heavy mesons in Rcfs. 
|1 5l — tl 9j | . In this work we go one step further and study the possible molecular system composed of a pair of P-wave 
■"si" ■ heavy mesons in the (0 + , 1 + ) doublet according to the classification of the heavy quark symmetry. 

The non-strange P-wave (0 + ,l + ) heavy mesons are very broad with a width around several hundred MeV [20j . 
Instead of forming a stable molecular state, the system composed of a pair of non-strange P-wave heavy meson decays 
' rapidly. Experimental identification of such a molecular state will be very difficult. The attractive interaction between 
■^j. . the meson pairs may lead to a possible threshold enhancement in the production cross section. 

In contrast, the P-wave charm-strange mesons D s o(2317) and -D s i(2460) lie below the DK and D*K threshold 
respectively. They are extremely narrow because their strong decays violate the isospin symmetry. The future 
■ experimental observation of the possible heavy molecular states composed of a pair of excited charm strange mesons 
may be feasible if they really exist. We study the charmonium-like system composed of a pair of excited charm strange 
mesons in this work. 

This paper is organized as follows. We review the formalism in Section [TT1 and present the results in Section ITnWl 
The last section is a short summary. 



II. FORMALISM 
Flavor wave functions 



We list the flavor wave functions of the possible molecular states composed of the P-wave (0 + , 1 + ) heavy doublet 
in Table Hill D* a denotes (Df, D* + , D*+) while D x denotes {D^,D+,D+ ). The neutral D* Q - D x system with the 
parameter c = ±1 corresponds to the positive and negative charge parity respectively. 
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TABLE I: The flavor wave function of the Dq — _D ( * and D± — D\ system. 



Do - DZ 



state wave function 

"i* WW 

D* -D* 
$° -^(D* °D* °-D* +D* - 
*° S ±(DfDf + D*+D* 



<K 
<h<> 

x s 



75*- n*+ 

d:+d*z 



D1-D1 



state 



wave function 



$**- 

$**° 

$**- 



d-di 
D+m 

DIM 
DtxDT 



TABLE II: The flavor wave function of the Dq 
negative charge parity respectively. 



Di system. The parameter c = ±1 for the Dq — Di system with positive and 



state 



$*+/$*+ 
$*-/$*- 



DS-Di 



wave function 



^{ir+m+cD+m ) 

^(d^Di+cdId;-) 

\[{D*°Dl + cD°Dq°) - (D*+D- + cD+D* ~)] 
+ cDlDt ) + {Dl+D- + cD+D* -)] 
-^(D'+DUcD+D*, ) 
^{D* Q °D: 1+ cDlD* M ) 
^(Dl+D^+cD+D*-) 

-L(d* + d: 1 + cd+d:c) 

^(DltDIi+cDfiD'to) 



B. Effective lagrangian 



With the help of the heavy quark symmetry and chiral symmetry, the strong interaction between the P-wave 
(0 + , 1 + ) heavy doublet reads 

C = ig'TrlSbj^A^S^+iP'TriSbV^V^- Pll ) ba S a ] 

+i\ / Tr[S b a^F^(p) ba S a ] + g' a Tr[S a aS a ], (1) 

where S represents the (0 + , 1 + ) doublet. Its matrix representation is 

S = 1(1 + ^7,75-^], (2) 
S = 7°^ t 7°- (3) 
At the leading order, the axial vector field reads 

Kb = l^d^ ~ ^) ab = j-d»V ab + . . . , (4) 



where 



71 + ^ - + K " 
K - R o _^ 
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p^ b and F fMU (p) a i > represent the vector meson field and its strength tensor 



a 
Pab 



l 9v yH 

V2 Qb ' 



V2 



(0"V - d v V») ab + ..., 



where g v = m p / with m p = 0.77 GeV and /„• = 0.132 GeV. V is the nonet vector meson matrices 



V = 



\ 



P 

K* 



K*° 




Similarly the scalar field a represents the scalar nonet. All the coupling constants g' , /3', A' and g' a are real. 
In our calculation we only need the effective lagrangian at the tree level 



Cd 1 d 1 v 



£>D* D*o 



V2g v l3'v^V^) ba D* b D* l 

-VWV'O^&aODu • D\ a ) + 2V2ig v X'(d IJ ,V v - W^D^D^, 
-y/2g v \'{d,y v - d u V^)bae a ^ f) vp{D lba Dll + D* b D\ aa ), 



t oylfiv 



C 



D 1 D 1 a 



Jit 

2g' a D* Qa Dlla, 
-2g' a {D la -D\ a )a. 



(6) 



None of the coupling constants g' ', A', g' a are known precisely although there exists some crude theoretical estimation 
[2ll ]. We allow the parameters involved in this work to vary around the values extracted from the QCD sum rule 
approach (QSR). 



C. Derivation of the effective potential 

We follow Refs. [TBI . [l6| to derive the effective potential of the heavy molecular system. Interested readers may 
consult Refs. [H], Il6l| for details. As usual, the monopole type form factor (FF) is introduced at every interaction 
vertex in order to account for the non-point-like structure effect of each interacting particle and cure the singularity 
of the effective potential. 

A 2 — w? 

= A^- (?) 

A is the phenomenological cutoff parameter. Generally A is expected to be larger than the exchanged meson mass 
and lies around 1-3 GeV. 

The effective potential in the coordinate space reads 

V(r) = 7i/ d q V(q)F(q) 2 e-^, (8) 

> Y(A,m,r), (9) 

► Z(A,m,r), (10) 



q 2 + m 2 

q 2 

q 2 + to 2 
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Y(A,m,r) = ^(e— - e"-) - ^e"-, (11) 

I Q Q 

Z(A,m,r) = - — — (r 2 — )Y(A,m,r) 

e -mr m 2 e -Ar£2 e ~Ar^ A £ -Ar A 2 



4^7" 47rr 87r 47rr 



with £ = VA 2 — m 2 . 

We collect the meson masses in Table IIIII 



TABLE III: The meson masses 



meson 


mass (GeV) 


mason 


mess (GeV) 


meson 


mass (GeV) 




2.4 




2.4 




2.317 


D? 


2.42 


Df 


2.42 




2.46 


P° 


0.77 


P ± 


0.77 




0.782 


4> 


1.020 


ir° 


0.135 




0.140 


V 


0.548 


a 


0.66 


/o(980) 


0.98 



III. THE Do - Do CASE 

In the Dq — Dq case, the pseudoscalar meson exchange is forbidden by parity and angular momentum conservation. 



$f and $s( < I ) °) states don't exist. For $ ± , <&°, the effective potential reads 

V{r)tj al = -\g 2 J' 2 [-Y(A : m p7 r)+Y(A, m ^r)}-g' 2 Y(A,m a ,r), (13) 

^(r)*L = -J.9, 2 /3' 2 [3V(A,m p ,r) + y(A,m LJ ,r)]-.g; 2 y(A,m a ,r), (14) 

^(r)?L = -^/^(A, m , r) - g' 2 Y(A, m /o ,r), (15) 



We use the MATSLISE package to solve Schrodinger equation with the effective potentials. We collect the variation 
of the binding energy E (in unit of MeV) and the root-mean-square radius r (in unit of fm) with the cutoff and the 
coupling constants in Tabic llVl 

As the coupling constants increase, the attraction becomes stronger. The cutoff parameter reflects the non-point- 
like structure of the interacting hadrons at each vertex. Its value is the hadronic size. In this work we assume the 
"reasonable" cutoff should be larger than the exchanged light meson mass and be around 1-3 GeV. 

Simply for comparison, we also collect the numerical results for the other possible molecular states in the same 
multiplct although their experimental observation may be difficult because of the broad width of the non-strange 
(0 + , 1 + ) charmed mesons. 

IV. THE Di - Di CASE 

The effective potential of the $** ± , $**°, $|*°, $** ± ! $**°, ®*f systems reads 

V(r)ToTai° lJ] = -\g 2 J 2 C{J)[-Y{A, m p , r) + Y(A, m u ,r)] - X' 2 g 2 v B( J)[-Z{A, m p , r) + Z(A, m u , r)] 

+ ^2-A(J)l-Z(A, m v , r)) + \z{A, m„, r))] - g' 2 C( J)Y(A, m a , r), (16) 

V(r)f2i ] = - \9ll3' 2 C(J) [3Y(A, m p ,r) + Y(A, m u , r)] - X' 2 g 2 B(J) [3Z(A, m p ,r) + Z(A, , r)] 

+( n A(J)[3Z(A, m v ,r)) + \z{A, m v ,r))] - g' 2 C(J)Y(A, m a , r), (17) 
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TABLE IV: The variation of the binding energy E (in unit of MeV) and the root-mean-square radius r rms (in unit of fm) with 
the cutoff and the coupling constants for the Dq — Dq system. 

P' = 0.84, g' a = 0.761 



states 


A 


E 


>rms 


$ 










1.6 


-8.1 


1.28 




1.7 


-14.4 


1.01 




1.8 


-21.9 


0.86 




1.9 


-30.3 


0.76 










P' = 0.98, g' a = 0.761 


states 


A 


E 




$ 








<E> 8 


1.4 


-9.2 


1.23 




1.45 


-14.0 


1.04 




1.5 


-19.6 


0.92 




1.55 


-25.8 


0.83 










P' = 1.12, g' a = 0.761 


states 


A 


E 


T*rms 


$ 








$8 


1.25 


-6.9 


1.39 




1.3 


-12.9 


1.09 




1.35 


-20.4 


0.92 




1.4 


-29.2 


0.81 




3.0 


-17.8 


0.87 




3.2 


-25.6 


0.75 




3.4 


-34.1 


0.66 



V(rp ota r = -f^A(J)Z(A,m v ,r), (18) 



nr) *r°/*rv] = _£_ A{J)z{Kmri ^ (19) 



^M?2 J1 = -\glP' 2 C{J)Y(k, m ^r)~2\i*glB{J)Z(k,m^r) 

+^A(J)Z(A, m v , r) - g' a 2 C(J)Y(A, m fo ,r), (20) 

6 J IT 

where A(J), B(J) and C(J) denote 

A(J) = Yl (lAi;lA 2 |J, ? n)(lA 3 ;lA 4 |J,m)l[ e t 1 -( ( fx6t 3 *)et 2 -(<rx e 'f*)], (21) 

Ai A2 A3 A4 

B{J) = (lAislA^mKlAaslA^J.m)!^ (22) 

Ai A2 A3 A4 

C(J) = J2 (lAi;lA 2 |J J m)(lA 3 ;lA4|J,m)[( e -t 1 -e1 3 *)(et 2 -ei 4 *)] J (23) 

Ai A2A3A4 

?b ?2, ?3, ?4 are the polarizations of the initial and final states, c.t.s denotes 

c.t.s = (if* • £Kef* • ^(ef 1 • ??) - (ef ■ <f)(^* ■ SX^ 2 ■ ef*) 

The values of -4(J), £>(</) and C(J) with different quantum numbers are listed in Table IVl 
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TABLE V: The values of A( J), B(J) and C(J) with different quantum numbers. 



A(J) B{J) C(J) 



We collect the variation of the binding energy E and the root-mean-square radius r with the cutoff and the coupling 
constants in Tables IVTll VIII With the pscudoscalar meson exchange force alone and g' — 0.80, there exists an isovector 
$** with J p = + , a $g* state with J p = 2 + , an isoscalar $** state with J p = 0+ and an isoscalar state with 
J p = 2 + . Increasing g' to 1.06, we can find the bound state solution for $** with J p = 1 + and $** with J p = 1 + 
besides the above mentioned bound states. With g' = 1.32, the above bound states still exist. We notice that 
corresponding cutoff A becomes smaller with the larger g' . 

Including all the exchange meson contributions, we list the numerical results in Table IVIII From Tables I VII and 
IVIII we note that the pseudoscalar meson exchange potential is dominant in the total effective potential. Thus, it is 
reasonable to consider pseudoscalar meson exchange potential only when studying whether there exists a bound state 
solution for the D\ — Di case. Meanwhile, for the B\ — B\ system, we list the results in Table IXll 



V. THE Do - D x CASE 



In the Dq — Di case, there are both direct and crossed scattering channels in the derivation of the effective potential 
in the momentum space. In the crossed channel, the mass difference go between the initial and final states (i.e., Dq 
and Di) should be kept. We introduce 



Mm = \Jm 2 - %, 
a = v/A 2 - gg, 

where the subscript "m" denotes the exchanged meson. Accordingly, 

a 2 — /r 



F(q) = 



A 2 — q 2 a 2 + q : 



2 



After the Fourier transformation 



-a 2 2 — >■ ~ Y (a,n,r), (24) 

q z — m z gg — q m z 

2 2 

— g = ~ 2 2 — * -Z(a,[i,r), (25) 

q z — m z gg — q m z 

the effective potential for the $* ± , $*°, $*°, $* ± , $*±, $*°, $*°, $*? systems read as 

^(0™aI /$ * ±,<] = -\g 2 J 2 [-Y{k, m p , r) + Y(A, m u ,r)) + c{ ^g 2 v X' 2 [-Z(a, Mp , r) + Z(a, ^,r)] 



V(r) 



*8 7*8 

Total 



^g 2 ,/3' 2 [3Y (A, m„ r) + y (A, to w , r)] + c{ ^g 2 A' 2 [3Z(a, fx p , r) + Z(a, ^ , r)] 

1 g' 2 1 ~i 

q J2 [3^(a,/x ff ,r) + -Z(a,^,r)]| - g£ F(A, r), 



(27) 



9 / 2 



Z(a,ft v ,r), 



(28) 



^ V J Total 



1/ 



■Z(a,ii v ,r), 



(29) 
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Virpl'T* 1 = -\g 2 J 2 Y(K, m+, r) + c{ p 2 ,X' 2 Z(a, H , r) - \j^Z{a, M „, r)} - g' 2 Y(A, m fo , r). (30) 

The Z?o — -Di system is very similar to the D — D* case and is particularly interesting since A (3872) is often 
speculated to be a D — D* molecular candidate. The only difference is that both components in the D$ — D\ system 
are extremely narrow P-wave states. We first focus on the pseudoscalar meson exchange, which is repulsive for the 
$*i state with negative charge parity. The J PC = 1 ++ state appears as shown in Tables IVIIIHXI By comparing 
the result listed in Table I Villi and that in Tabic IIX1 one notices that the pseudoscalar meson exchange is dominant 
in the Dq — D\ system, which shows that it is reasonable to consider the pseudoscalar meson exchange potential only 
when we investigate whether there exists the bound state solution for the Dq — D\ system. The result for the — B\ 
system corresponds to the pseudoscalar meson exchange only. 

VI. CONCLUSION 

Both Z? s o(2317) and -D s i(2460) lie below the DK and D*K threshold respectively. They are extremely narrow. 
The possible molecular states composed of the £> s o(2317) and Z? s i(2460) may be observable experimentally if they 
really exist. In this work we have studied such systems carefully. As a byproduct, we collect the numerical results for 
the bottomonium-like analogue in the appendix. 

One should be cautious that our numerical results are quite sensitive to the values of the hadronic coupling constants. 
The values are larger than (or around the upper bound of) those derived from the crude estimate with the light cone 
QCD sum rule approach [2l| . Future lattice QCD simulations may help extract these coupling constants more precisely. 
Since the hadronic coupling constants are not known well, we allow them to vary. As shown in the numerical result, 
the binding energy is also sensitive to the value of the cutoff introduced in the form factor. Thus, further study and 
improvement of the potential model are still desirable. 

Here, we need to emphasize that a monopole form factor is introduced in the numerical calculation of this work. 
In fact, there are many types of form factor, such as the dipole form factor. When taking the other type of the form 
factor, the qualitative conclusion keeps the same as that obtained in this work. Both the form factor and the cutoff 
are necessary and important for the hadronic system since the components are not point-like particles. They are 
hadrons with internal structure. When dealing with the loosely bound heavy molecular states, only the relatively soft 
degree of freedom is expected to play the dominant role. The exchanged soft mesons should not "see" the quark/gluon 
structure of the heavy meson. That's the physical meaning of the form factor and the cutoff. 

So long as these couplings are big enough, there may even appear deeply bound states including radial and orbital 
excitations. However they are no more the "conventional" molecular states, which are loosely bound with a typical 
binding energy around several to several tens MeV and a radius around 1.5-3 fm. Therefore we do not list numerical 
results for the deeply bound cases in this work. 

From our calculation there may exist two loosely bound ++ charmonium-likc states, the first of which is composed 
of the D s oD s q pair and lies around 4.61 GeV. The other one is around 4.9 GeV and composed of the D s \ and D s \ 
pair. There exists the 2 ++ D s iD s i state, which lies around 4.9 GeV. The 1 ++ state around 4.75 GeV is composed of 
the D s q and D s \ pair. This state is very interesting because of its similarity to A (3872). 

The dominant decay modes of the above states are the open-charm modes D), ' D\ '. The other characteristic 
decay modes are the hidden-charm modes J/il)<j>, r] c r]', ?7c/o(980), Xcjv\ Xcj/o(980), tp'cf), ip"cj), rj c (2S)r]' , etc. for the 
possible C = + molecular states. One may easily exhaust the possible final states according to the C/P parity and 
angular momentum conservation and kincmatical considerations. These states may be significantly narrower than 
the conventional charmonium around the same mass region because of their molecular nature. However, their widths 
should be larger than those of A (3872) due to much larger phase space and more decay modes. 

These states might be produced from B or B s decays if kinematically allowed. Those states with J PC = ++ , 2 ++ 
may be produced from the two photon fusion process at the e + e~ collider at B factories. The other possible facilities 
to look for them are RHIC, Tevatron and LHCb. Investigations of these states may help us understand the puzzling 
A (3872) state. 
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TABLE VI: The variation of the binding energy E (in unit of MeV) and the root-mean-square radius r rms (in unit of fm) with 
the cutoff and the coupling constant for the D\ — D\ system when only the pseudoscalar meson exchange is considered. Here, 
we scan the cutoff range A < 3.2 GeV. 

g' = 0-80, P' = 0, A' = 0, g ' a = 



state 




J p = l+ 


J p = 2+ 


A E Trms 


A E Trms 


A E T rrns 


$** 


1.9 -5.4 1.33 

2.0 -10.6 0.98 

2.1 -17.9 0.78 

2.2 -27.9 0.65 


- 


- 


$** 






1.1 -4.5 1.49 

1.2 -13.6 0.93 

1.3 -28.9 0.69 




2.7 -6.58 1.13 

2.8 -18.5 0.70 

2.9 -36.4 0.52 












2.7 -7.8 1.03 

2.8 -20.5 0.66 

2.9 -39.2 0.50 



1.06, P' = 0, A' = 0, 



state 


J p = 0+ 


J p = l + 


J p = 2+ 


A E T r m S 


A E T r ms 


A E T r ms 


$** 


1.1 -4.0 1.58 

1.2 -9.9 1.08 

1.3 -19.0 0.82 

1.4 -31.9 0.67 


2.1 -3.4 1.62 

2.2 -7.4 1.14 

2.3 -13.3 0.88 

2.4 -21.3 0.72 




$8* 






0.75 -3.6 1.72 
0.8 -8.0 1.23 
0.9 -23.6 0.81 




1.9 -4.3 1.39 

2.0 -18.9 0.71 

2.1 -43.4 0.50 


3.0 -11.7 0.86 

3.1 -25.6 0.60 

3.2 -44.8 0.47 










1.9 -5.1 1.28 

2.0 -20.5 0.68 

2.1 -45.7 0.49 


3 ' = 1.32, /?' = 0, A' = 0, g' a = 


state 


J p = 0+ 


J^ = l+ 


J p = 2+ 


A £/ Trins 


A E T rms 


A E Trms 


$" 


0.75 -3.0 1.86 
0.85 -10.3 1.11 
0.95 -22.8 0.82 
1.05 -41.3 0.65 


1.35 -2.7 1.85 
1.50 -10.0 1.04 
1.65 -23.1 0.73 
1.75 -35.9 0.61 










0.65 -12.0 1.11 
0.70 -22.0 0.88 
0.75 -35.9 0.73 




1.55 -8.4 1.03 
1.60 -18.5 0.73 
1.65 -32.4 0.58 


2.25 -9.7 0.95 
2.30 -16.9 0.74 
2.40 -37.1 0.53 










1.50 -2.6 1.78 
1.55 -9.3 0.98 
1.60 -19.8 0.71 
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TABLE VII: The variation of the binding energy E (in unit of MeV) and the root-mean-square radius r rms (in unit of fm) 
with the cutoff and the coupling constants for the Di — Di system. Here, we scan the cutoff range A < 3.2 GeV. 

g' = 0.80, 13' = 0.84, A' = 0.42, g' a = 0.761 



state 


J p = 0+ 


J^ = l+ 


J p =2+ 


A E Trms 


A E Trms 


A E frms 


$** 


1.9 -7.1 1.18 

2.0 -13.0 0.91 

2.1 -21.0 0.74 

2.2 -31.7 0.62 


- 




*r 






1.0 -5.7 1.41 
1.05 -11.6 1.06 

1.1 -19.5 0.87 




2.7 -6.58 1.13 

2.8 -18.5 0.70 

2.9 -36.4 0.52 












2.6 -3.2 1.70 
2.9 -10.0 1.03 
3.2 -18.8 0.79 
3.5 -28.6 0.66 



g = 1.06, P' = 0.98, A' = 0.49, g' a = 0.761 





J p = 0+ 


J^ = l+ 


J p =2+ 


state 


A 


E 


' rvns 


A 


E 


frms 


A 


E 


frms 


$** 


1.1 


-4.1 


1.57 


2.1 


-5.1 


1.35 










1.2 


-10.1 


1.07 


2.2 


-9.9 


1.01 










1.3 


-19.3 


0.82 


2.3 


-16.6 


0.81 










1.4 


-32.4 


0.66 


2.4 


-25.5 


0.67 






















0.8 


-8.4 


1.21 
















0.85 


-16.9 


0.93 
















0.9 


-29.7 


0.76 




1.9 


-4.3 


1.39 


2.9 


-3.1 


1.61 










1.95 


-10.4 


0.92 


3.0 


-11.7 


0.86 










2.0 


-18.9 


0.71 


3.1 


-25.6 


0.60 










2.05 


-29.9 


0.58 


3.15 


-34.5 


0.53 






















1.75 


-3.6 


1.60 
















1.8 


-6.9 


1.20 
















1.9 


-16.1 


0.83 
















2.0 


-28.1 


0.66 



g' = 1.32, /3' = 1.12, A' = 0.56, g' a = 0.761 



state 


J p = 0+ 


J^ = l+ 


J p =2+ 


A E Trms 


A E Trms 


A E Trms 


$** 


0.8 -6.1 1.38 
0.85 -10.3 1.12 
0.9 -15.7 0.95 
0.95 -22.5 0.82 


1.4 -4.9 1.42 

1.5 -10.6 1.02 

1.6 -19.0 0.80 

1.7 -30.5 0.66 




®r 






0.8 -55.1 0.62 
0.9 -118.0 0.48 
1.0 -213.9 0.39 




1.55 -8.4 1.03 
1.60 -18.7 0.73 
1.65 -32.4 0.58 


2.2 -4.4 1.42 
2.25 -9.7 0.95 

2.3 -16.9 0.74 










1.45 -5.4 1.32 
1.50 -12.5 0.92 
1.55 -22.0 0.73 
1.60 -33.7 0.62 
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TABLE VIII: The variation of the binding energy E (in unit of MeV) and the root-mean-square radius r rms (in unit of fm) 
with the cutoff and the coupling constant for the Dq — D\ system when only the pseudoscalar meson exchange is considered. 



g' = 0.80, 0' = 0, A' = 0, g' a = 




c = +l 


c= -1 


states 


A E Trms 


A E Trrns 








$8 


1.1 -4.4 1.51 

1.2 -13.4 0.94 

1.3 -28.6 0.69 

1.4 -51.3 0.55 






2.7 -5.2 1.27 

2.8 -16.0 0.75 
2.85 -23.6 0.63 

2.9 -32.6 0.55 





g' = 1.06,(3' = 0, \' = 0,9^ = 





c= +1 


c = -1 


states 


A £/ T^rms 


A E Trms 


$* 




2.2 -7.1 1.17 

2.3 -12.8 0.90 

2.4 -20.7 0.73 

2.5 -30.9 0.61 




0.8 -8.0 1.24 
0.85 -14.5 0.98 
0.9 -23.5 0.81 
0.95 -35.3 0.69 








3 -9.0 0.98 
3.05 -14.5 0.78 
3.1 -21.4 0.66 
3.15 -29.6 0.57 




1.9 -4.2 1.43 
1.95 -10.1 0.95 
2.0 -18.3 0.73 
2.05 -29.0 0.60 




g' = 1.32, P' = 0, A' = 0, = 




c=+l 


c = -1 


states 


A E Trms 


A E Trms 






1.5 -9.7 1.06 

1.6 -17.7 0.82 

1.7 -28.6 0.67 




0.8 -54.2 0.63 
0.825 -65.4 0.58 
0.85 -77.9 0.55 
0.875 -92.0 0.51 








2.25 -7.9 1.05 

2.3 -14.5 0.80 
2.35 -23.0 0.65 

2.4 -33.5 0.55 




1.55 -8.9 1.02 

1.6 -19.1 0.73 
1.65 -32.9 0.59 

1.7 -50.5 0.49 





DPFIHE under Grants No. 20090211120029, NCET under Grants No. NCET-10-0442). 
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TABLE IX: The variation of the binding energy E (in unit of MeV) and the root-mean-square radius r rms (in unit of fm) with 
the cutoff and the coupling constants for the Dq — D\ system. 



9 — in- 


SO, p 


= 0.84, 


A' = 






c = +l 


c= -1 


states 


A 


E 


Trms 


A E Trrns 


$* 






rfi* 
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-5.6 


1.42 






1.05 


-11.5 


1.07 






1.1 
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0.88 






1.15 
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^ si 


2.7 


-4.0 
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3.0 
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3.3 
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0.79 
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-28.3 


0.67 
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Appendix A: Possible molecular states composed of a pair of excited bottom-strange mesons 

We collect the numerical results for the bottomonium-like system composed of a pair of excited bottom-strange 
mesons in the appendix. Since neither B s q nor B s i is observed experimentally, we follow Ref. p2| and use the mass 
values m Bo (J P = 0+) = 5.627 GeV, m Bl (J p = 1+) = 5.674 GeV, m Bs0 (J p = 0+) = 5.718 GcV and m Bsl (J p = 
1+) = 5.765 GeV. 

TABLE X: The variation of the binding energy E (in unit of MeV) and the root-mean-square radius r rms (in unit of fm) with 
the cutoff and the coupling constants for the Bg — Bq system. 



P' = 0.84, g' a = 0.761 



states 


A 


E 


Trms 


$ 








$8 


1.15 


-3.0 


1.41 




1.2 


-6.5 


1.04 




1.25 


-11.1 


0.85 




1.35 


-23.4 


0.66 




2.2 


-5.4 


1.00 




2.3 


-8.3 


0.84 




2.4 


-11.6 


0.73 


/?' = 0.98, g' a = 0.761 


states 


A 


E 














1.11 


-6.0 


1.08 




1.15 


-10.8 


0.87 




1.19 


-16.8 


0.75 




1.23 


-23.9 


0.66 




1.9 


-5.9 


0.98 




2.0 


-10.5 


0.78 




2.1 


-16.2 


0.66 




2.2 


-22.8 


0.58 


= 1.12, g' a =0.761 


states 


A 


E 


Trms 


$ ±,u 










1.05 


-5.6 


1.12 




1.1 


-13.3 


0.82 




1.15 


-24.0 


0.67 




1.2 


-37.3 


0.58 




1.7 


-5.0 


1.06 




1.75 


-7.9 


0.88 




1.8 


-11.3 


0.77 




1.85 


-15.2 


0.69 
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TABLE XI: The variation of the binding energy E (in unit of MeV) and the root-mean-square radius r rms (in unit of fm) with 
the cutoff and the coupling constant for the B\ — B\ system when only the pseudoscalar meson exchange is considered. The 
Bi — Bi system is easier to form a bound state than the D\ — D± case. In this table, we only give the result for the B\ — B\ 
system with the typical coupling constant g' = 0.80, 1.06. We scan the cutoff range A < 3.1 GeV. 







/ 

9 


= 0.80, P' 


= 0, A' = 


0, g'a 


= 








J p = + 


J p = l+ 




state 


A 


E 


frnns 


A 


E 


frms 


A 


E 


V*rms 


$** 


1.0 


-6.4 


0.94 


1.6 


-1.7 


1.52 










1.1 


-12.2 


0.73 


1.8 


-6.8 


0.84 


- 


- 


- 




1.2 


-20.3 


0.60 


2.0 


-16.5 


0.58 










1.3 


-31.2 


0.50 


2.2 


-32.2 


0.44 






















0.8 
0.85 


-13.9 
-20.8 


0.74 
0.64 
















0.9 


-29.7 


0.56 




1.7 


-6.2 


0.81 


2.45 


-3.4 


1.03 










1.8 


-17.4 


0.53 


2.55 


-9.3 


0.65 










1.9 


-34.5 


0.40 


2.65 


-18.2 


0.49 








<E>*i* 

^ si 














1.65 
1.7 
1.8 


-3.0 
-6.7 
-18.3 


1.10 
0.78 
0.51 






9 


= 1.06, 0' 


= 0, A' = 


0, g'a 


= 








J p = 0+ 


J p = 1+ 




state 


A 


E 


Trms 


A 


E 


Trms 


A 


E 


Trms 


$** 


0.8 


-18.2 


0.69 


1.1 


-6.1 


0.95 










0.85 


-24.8 


0.61 


1.2 


-11.2 


0.74 










0.9 


-32.8 


0.55 


1.3 


-18.4 


0.61 






















0.8 
0.825 


-70.2 
-81.8 


0.44 
0.41 
















0.85 


-94.7 


0.39 


^ s 


1.3 


-3.8 


1.02 


1.8 


-5.0 


0.88 










1.35 


-9.5 


0.70 


1.9 


-14.6 


0.56 










1.4 


-17.6 


0.55 


2.0 


-29.3 


0.42 






















1.3 
1.35 
1.4 


-4.1 
-10.0 
-18.3 


0.98 
0.68 
0.54 
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TABLE XII: The variation of the binding energy E (in unit of MeV) and the root-mean-square radius r rms (in unit of fm) 
with the cutoff and the coupling constant for the Bq — B\ system when only the pseudoscalar meson exchange is considered. 

g = 0.80, /?' = 0, A' = 0, g' tT = 





c=+l 


c = -1 


states 


A Trms 


A E Trms 


$* 




1.6 -1.8 1.50 
1.8 -6.9 0.84 
2.0 -16.5 0.58 
2.2 -32.2 0.44 


$* 


0.8 -14.3 0.74 
0.825 -17.6 0.68 
0.875 -25.6 0.59 

0.9 -30.3 0.56 




^ s 




2.45 -3.5 1.02 
2.55 -9.5 0.65 
2.65 -18.4 0.49 
2.75 -30.4 0.40 


^ si 


1.65 -3.0 1.12 
1.75 -11.7 0.62 
1.8 -18.1 0.52 
1.85 -26.0 0.45 




g' = 1.06, 13' = 0, A' = 0, g' a = 




c=+l 


c = -1 


states 


A E T'rms 


A E Trms 


$* 




1.05 -4.4 1.09 
1.25 -14.8 0.66 
1.35 -23.2 0.56 
1.45 -34.2 0.48 




0.8 -71.4 0.44 
0.825 -83.0 0.41 
0.85 -95.9 0.39 
0.875 -110.1 0.37 




<&: 




1.75 -2.3 1.26 
1.85 -9.6 0.67 
1.95 -21.8 0.48 
2.05 -39.3 0.38 




1.3 -4.1 0.98 
1.35 -9.9 0.68 

1.4 -18.2 0.54 
1.45 -29.1 0.45 




g' = 1.32, p = 0, A' = 0, = 




c= +1 


c = -1 


states 


A E Trms 


A E Trms 


$* 




0.8 -6.9 0.96 
0.9 -14.2 0.73 

1.0 -24.9 0.58 

1.1 -39.5 0.49 




0.8 -177.2 0.33 
0.81 -187.1 0.32 
0.82 -197.3 0.31 
0.83 -208.1 0.31 








1.4 -4.6 0.94 

1.5 -9.7 0.68 

1.6 -25.8 0.46 
1.65 -36.8 0.40 




1.1 -3.1 1.13 
1.15 -10.0 0.70 

1.2 -20.9 0.53 
1.25 -35.8 0.43 





